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The aim of this paper is two-fold: in probing the statistical mechanical properties of 
interacting quantum fields, and in providing a field theoretical justification for a stochastic 
source term in the Boltzmann equation. We start with the formulation of quantum field 
theory in terms of the Schwinger - Dyson equations for the correlation functions, which we 
describe by a closed-time-path master (n = ooPI) effective action. When the hierarchy is 
truncated, one obtains the ordinary closed-system of correlation functions up to a certain 
order, and from the nPI effective action, a set of time-reversal invariant equations of motion. 
But when the effect of the higher order correlation functions is included (through e.g., causal 
factorization- molecular chaos - conditions, which we call 'slaving'), in the form of a correla- 
tion noise, the dynamics of the lower order correlations shows dissipative features, as familiar 
in the field-theory version of Boltzmann equation. We show that fluctuation-dissipation re- 
lations exist for such effectively open systems, and use them to show that such a stochastic 
term, which explicitly introduces quantum fluctuations on the lower order correlation func- 
tions, necessarily accompanies the dissipative term, thus leading to a Boltzmann-Langevin 
equation which depicts both the dissipative and stochastic dynamics of correlation functions 
in quantum field theory. 

I. INTRODUCTION 

The main result of this paper is a derivation of the Boltzmann - Langevin equation as the correct description 
of the kinetic limit of quantum field theory. We use for this derivation a generalized formulation of quantum 
field theory in terms of the Schwinger - Dyson equations but supplemented with stochastic terms which 
explicitly introduce the effects of quantum fluctuations on low order correlation functions. 

The significance of such an inquiry is two-fold: In probing the statistical mechanical properties of inter- 
acting quantum fields, and in providing a field theoretical justification for a stochastic source term in the 
Boltzmann equation. The former has been investigated mainly for quantum mechanical, not field-theoretical, 
systems |j],|2| (see however, and the latter primarily for a classical gas Extending previous studies 

to quantum fields is essential in the establishment of a quantum field theory of nonequilibrium processes. 
Previous work on this subject showed how the Boltzmann equation can be derived from first-principles 

in quantum field theory and investigated its dissipative properties. In this paper we focus on the fluctuation 
and noise aspects in the derivation of a stochastic Boltzmann equation from quantum field theory [jl5|— ]5~7[] . 

Fluctuations in Composite Operators 



There are a variety of problems in nonequilibrium field theory which are most naturally described in terms 
of the unfolding of composite operators, the familiar lowest order ones being the particle number and the 
energy momentum densities and their fluctuations. The usual approach to these problems assumes that these 
operators have small fluctuations around their expectation values, in which case they can be expressed in 
terms of the Green functions of the theory. However, when fluctuations are large, typically when correlations 
among several particles are important, this approximation breaks down. 

A familiar example is critical phenomena: by choosing a suitable order parameter to describe the different 
phases, one can obtain a wealth of information on the phase diagram of a system. But to study the dynamics 
of a phase transition, especially in the regime where fluctuations get large, the single order parameter must 
be replaced by a locally defined field obeying a stochastic equation of motion, for example, a time - dependent 
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Ginzburg - Landau equation with noise (which though often put in by hand, should in theory be derived 
from fluctuation - dissipation considerations). The same phenomenon occurs more generally in effective 
field theories, where the light fields are randomized by the back reaction from the heavy fields pl| , and in 
semiclassical theories, where the classical field (for example, the gravitational field in the Early Universe) is 
subject to random driving forces from activities in the quantum field, such as particle creation |js[ ]. 

The influence of noise on the classical dynamics of a quantum system is discussed at length by Gell-Mann 
and Hartlc [ p0[ | ; the conversion of quantum fluctuations to classical noise is discussed by a number of authors 
pl| p5| . This scheme was also used b y u s for the study of decoherence of correlation histories and correlation 
noise in interacting field theories @^^] and related work by others using projection operators P^| . 

Therefore in the more general context illustrated by these examples, we shall explore models where the 
dynamics is described not in terms of the actual Green functions G ab (we use closed - time - path (CTP) 
techniques and notation, described more fully in the Appendix ^9|) but a stochastic kernel G ofc , whose 
expectation value reduces to G ab , while its fluctuations reproduce the quantum fluctuations in the binary 
products of field operators. 

More concretely, consider a theory of a scalar field <fi a (we use a condensed notation where the index a 
denotes both a space - time point and one or the other branch of the time path - see Appendix). The CTP 
action is S — S [(f) 1 ] — S* [</> 2 ] . Introduce the generating functional 
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This suggests viewing the stochastic kernel G ab as a Gaussian process defined (formally) by the relationships 

(G ab ) = {<p a <p b ) ■ {G ab G cd ) = (</>>V c d ) (1-4) 

Or else, calling 

G ab = G ab + A ab 



(A ab ) = 0; (A ab A cd ) 



-4/ 
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(1.6) 



K=0 



To turn the intuitive ansatz Eqs. (1.4) and (1.6) into a rigorous formalism we must deal with the obvious 
fact that we are manipulating complex expressions; in particular, it is not clear the As define a stochastic 
process at all. However, for our present purposes it will prove enough to deal with the propagators as if they 
were real quantities. The reason is that we are primarily concerned with the large occupation numbers or 
semiclassical limit, where the propagators do become real. We will see that this prescription will be sufficient 
to extract unambiguous results from the formal manipulations below. 



Stochastic Boltzmann Equation 

We can define the process A ab also in terms of a stochastic equation of motion, 
transform of W, the so - called 2 particle irreducible (2PI) effective action (EA) 



r [G ab ] = W [K* ab ] 
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We have the identities 



ST 
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0: 
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(I. 



the first of which is just the Schwinger 
following equations of motion for G ab 



Dyson equation for the propagators; we therefore propose the 
ST -1 



5G ab 2 

where K a b is a stochastic nonlocal Gaussian source defined by 



(K a b) = 0; 



(KabKcd) = 4i 



<5 2 r 



1 1 



6G ab SG cd 



(1.9) 



(1.10) 



If we linearize Eq. Jl.9|) around G, then the correlation Eq. ( 1. 10 ) for k implies Eq. (1.6) for A. Consistent 
with our recipe of handling G as if it wer e re al we should treat k also as if it were a real source. 

It is well known that the noiseless Eq. (L9) can be used as a basis for the derivation of transport equations 
in the near equilibrium limit. Indeed, for a type theory, the resulting equation is simply the Boltzmann 
equation for a distribution function / defined from the Wigner transform of G ab (details arc given below). 
We shall show in this paper that the full stochastic equation (L9) leads, in the same limit, to a Boltzmann - 
Langevin equation, thus providing the microscopic basis for this equation in a manifestly relativistic quantum 
field theory. 

Let us first examine some consequences of Eq. ([.10). For a free field theory, we can compute the 2PI EA 
explicitly (derivation in Section V) 



r [G ab ] = -1 In [DetG] - l -c ah (-□ + m 2 ) G ab (x, x) 
where c a b is the CTP metric tensor (see the Appendix). We immediately find 

s 2 r 



6G ab 6G cd 



therefore 



(A ab A cd ) 



<5 2 r 



G ac G ti 



_SG ab SG cd . 

an eminently sensible result. Observe that the stochastic source does not vanish in this case, rather 

(nabK c d) = G a ^G di + G d aG b ^ 

However 



-iCac (— □ 



(1.11) 



(1.12) 



(1.13) 



(1.14) 



(1.15) 



does vanish on mass - shell. Therefore, when we take the kinetic theory limit, we shall find that for a free 
theory, there are no on - shell fluctuations of the distribution function. For an interacting theory this is no 
longer the case. 

The physical reason for this different behavior is that the evolution of the distribution function for an 
interacting theory is dissipative, and therefore basic statistical mechanics considerations call for the presence 
of fluctuations |3(|. Indeed it is this kind of consideration which led us to think about a Boltzmann - 
Langevin equation in the first place. This is fine if one takes a statistical mechanical viewpoint, but one 
is used to the idea that quantum field theories are unitary and complete with no information loss, so how 
could one see dissipation or noise? 



3 



In field theory there is a particular derivation of the self consistent dynamics for Green functions which 
resolves this puzzle, namely when the Dyson equations are derived from the variation of a nonlocal action 
functional, the two-particle irreducible effective action (2PI-EA). This was originally introduced |3l[] as a 
convenient way to perform nonperturbative resummation of several Feynman graphs. When cast in the 
Schwinger - Keldysh "closed time path" (CTP) formulation p9[ , it guarantees real and causal evolution 
equations for the Green functions of the theory. It is conceptually clear if one begins with a " master" effective 
action (MEA) j27j] where all Green functions of the theory appear as arguments, and then systematically 
eliminate all higher-than-two point functions to arrive at the 2PIEA. 

To us, the correct approach is to view the two point functions as an open system |32]-|34| truncated from 
the hierarchy of correlation functions obeying the set of Schwinger-Dyson equations but interacting with 
an environment of higher irreducible correlations, whose averaged effect brings about dissipation and whose 
fluctuations give rise to the correlation noise p7y . This is the conceptual basis of our program. 

Fluctuations and Dissipation 



It has long been known in statistical physics that the equilibrium state is far from being static; quite the 
opposite, it is the fluctuations around equilibrium which underlie and give meaning to such phenomena as 
Brownian motion p3[ and transport processes, and determine the responses (such as the heat capacity and 
susceptibility) of the system in equilibrium. The condition that equilibrium constantly reproduces itself in 
the course of all these activities means that the equilibrium state is closely related both to the structure of the 
equilibrium fluctuations and to the dynamical processes by which equilibrium sustains itself; these simple but 
deep relations are embodied in the so-called fluctuation - dissipation theorems: If a fluctuating system is to 
persist in the neighborhood of a given equilibrium state, then the overall dissipative processes in the system 
(due mainly to its interaction with its environment) are determined. Vice versa, if the dissipative processes 
are known, then we may describe the properties of equilibrium fluctuations without detailed knowledge of 
the system's microscopic structure. This is the aspect of the fluctuation - dissipation relations which guided 
Einstein in his pioneering analysis of the corpuscular structure of matter |35| ] , Nyquist in his stochastic theory 
of electric resistivity , and Landau and Lifshitz to the theory of hydrodynamical fluctuations p7| . 

These ideas apply to systems described by a few macroscopic variables, as well as to systems described 
by an infinite number of degrees of freedom such as a few long wavelength modes (as in hydrodynamics) or 
a single particle distribution function (as in kinetic theory). In this later case, the dynamics is described by 
a dissipative Boltzmann equation, and thereby we are to expect that there will be nontrivial fluctuations in 
equilibrium. The stochastic properties of the Boltzmann equation has been discussed by Zwanzig, Kac and 
Logan, and others P,p8[. 

The Boltzmann equation can be retrieved also as a description of the long range, near equilibrium dynamics 
of field theories JlO| . In this kinetic theory regime, where there is a clear separation of microscopic and 
macroscopic scales the field may be described in terms of quasi-particles, whose distribution function obeys 
a Boltzmann equation Jl^,^]. Formally, the one particle distribution function is introduced as a partial 
Fourier transform of a suitable Green function of the field. The same arguments which lead to a fluctuating 
Boltzmann equation in the general case, lead us to expect fluctuations in this limiting case of field theory as 
well. 

The end result of our investigation is a highly nonlinear, explicitly stochastic Dyson equation for the Green 
functions. By going to the kinetic theory limit, we derive a stochastic Boltzmann equation, and the resulting 
noise may be compared with that required by the fluctuation - dissipation relation. Here we see clearly this 
contrast between the predictions of field theory with and without statistical physics considerations. 

In this paper, we shall concentrate on the issue of what kind of fluctuations may be convincingly derived 
from the 2PI-CTP-EA for Green functions, and how they compare to the fluctuation - dissipation noise in 
the kinetic theory limit. Given the complexity of the subject, we shall adopt a line of development which 
favors at least in the beginning ease of understanding over completeness. That is, instead of starting from 
the master effective action of n point functions and work our way down in a systematic way, we shall begin 
with the Boltzmann equation for one -particle distributions and work our way up. 

In the next section, we present briefly the fluctuation - dissipation theorem in a nonrelativistic context, 
and use it to derive the fluctuating Boltzmann equation. The discussion, kept at the classical level, simply 
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reviews well established results in the theory of the Boltzmann equation. Section III reviews the basic tenets 
of nonequilibrium quantum field theory as it concerns the dynamics of correlations, and the retrieval of the 
Boltzmann equation therefrom. We refrain from using functional methods, so as to keep the discussion as 
intuitive as possible. 

Section IV discusses how the functional derivation of the Schwinger - Dyson hierarchy suggests that these 
equations ought to be enlarged to include stochastic terms. By going through the kinetic limit we use these 
results to establish a comparison with the purely classical results of Section II. 

Our investigation into the physical origin of noise and dissipation in the dynamics of the two point func- 
tions shows that in the final analysis this is an effective dynamics, obtained from averaging out the higher 
correlations. This point is made most explicit in the approach whereby the 2PI EA for the correlations is 
obtained through truncation of the master effective action, this being the formal functional whose variations 
generate the full Schwinger - Dyson hierarchy. In Section V, we briefly discuss the definition and construction 
of the master effective action, the relationship of truncation to common approximation schemes, and present 
explicitly the calculation leading to the dynamics of the two point functions at three loops accuracy |^7],[l0| ■ 

In the last section we give a brief discussion of the meaning of our results and possible implications on 
renormalization group theory. 



II. STOCHASTIC BOLTZMANN EQUATION FROM FDT 

As a primer, we wish to introduce the fluctuation-dissipation theorem (FDT) or relation (FDR) in its 
simplest yet complete form, and apply it to derive the stochastic Boltzmann equation so as to clarify its 
physical content. There are many different versions p9| : It could be taken to mean the formulae relating 
dissipative coefficients to time integrals of correlation functions (sometimes called the "Landau - Lifshitz 
FDT) or the relations between the susceptibility and the space integral of the correlation function. In this 
paper, the fluctuation-dissipation theorem addresses the relation between the dissipative coefficients of the 
effective open system and the auto- correlation of random forces acting on the system, as illustrated below. 



A. Fluctuation - dissipation theorem (FDT) 

The simplest setting for the FDT is a homogeneous system described by variables x l . The thermody- 
namics is encoded in the form of the entropy S (x l ). The thermodynamic fluxes are the derivatives x l , and 
the thermodynamic forces are the components of the gradient of the entropy 

*-* < IL1 » 

The dynamics is given by 

x^-fiFj+f (11.2) 

The first term describes the mean regression of the system towards a local entropy maximum, j lJ being the 
dissipative coefficient or function, and the second term describes the random microscopic fluctuations induced 
by its interaction with an environment. Near equilibrium, we also have the phenomenological relations for 
linear response 

Fi = djiP (II.3) 

where Cjj is a nonsingular matrix. 

In a classical theory, the equal time statistics of fluctuations is determined by Einstein's law 

{x l {t)F 3 (t))=5) (II.4) 

Take a derivative to find 
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= c jk {({-l U F l+ f) x k ) + (x* (- 7 H F; +/)>} 



(II.5) 



If the noise is gaussian, 



and white 



then 



(*' (t)j k (*)) 



dt 



, 8x l (t) 



Sj l (f ) 



(t)i*(*)> = 2" tt 



(j l (*')/(*)) 



(II.6) 



(II.7) 



From Eq. ([1. 5) and (II. 4) we find the noise-noise auto-correlation function v %k is related to the symmetrized 
dissipative function 7 l/c by 



(II. 



which is the FDT in a simple classical formulation^]. 

In the case of a one - dimensional system, the above argument can be simplified even further because there 
is only one variable x, and 7, c, v are simply constants. In equilibrium, we have (a; 2 ) = c^ 1 . On the other 
hand, the late time solution of the equations of motion reads 



x(t)= J du e^ c(t ' u) j(u) 
which implies (a; 2 ) = is/2"/c. Thus v = 27, in agreement with Eq. (|X 



B. Boltzmann equation for a classical relativistic gas 

We shall apply the theory above to a dilute gas of relativistic classical particles |Q. The system is 
described by its one particle distribution function /(X, A), where X is a position variable, and k is a 
momentum variable. Momentum is assumed to lie on a mass shell k 2 + M 2 — 0. (We use the MTW 
convention, with signature -H — I — h for the background metric f45|| ) and have positive energy fc° > 0. In other 
words, given a spatial element ci£ M = n^dT} and a momentum space element d k, the number of particles 
with momentum k lying within that phase space volume element is 

dn = -4tt/ {X, k) 9 (fc°) 6 (k 2 + M 2 ) k^n^ dS-^j (II.9) 

(2tt) 

The dynamics of the distribution function is given by the Boltzmann equation, which we give in a notation 
adapted to our later needs, and for the time being without the sought-after stochastic terms 

k^-^-f(k)=I col (X,k) (11.10) 



: To be concrete, this is the FDT of the second kind in the classification of Ref. Ell. The FDT of the first kind is 
further discussed in Ref. |l2|]. Also observe that we are only concerned with small deviations from equilibrium; FDT's 
valid arbitrarily far from equilibrium are discussed in Ref. [faja. 
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(11.11) 



! = {[! + / (ps)] [1 + / (*)] f (p x ) / (pa) - [1 + / (p x )] [1 + / (pa)] / (p 3 ) / (fc)} 
The entropy flux is given by 
d 4 p 



S>* (X) = 4tt 



(2tt)' 



-0 (p°) 5 (p 2 + M 2 ) p" {[1 + / (p)] In [1 + / (p)] - / (p) In / (p)} 



(11.12) 



(11.13) 



while the entropy itself 5 1 is (minus) the integral of the flux over a Cauchy surface. Now consider a small 
deviation from the equilibrium distribution 



/ = fea + Sf 



(11.14) 



i 



bPp - 1 

corresponding to the same particle and energy fluxes 

d 4 p 



(2*)< 



p°) S (p 2 + M 2 ) p^«5/ (p) = 



(11.15) 



(11.16) 



d 4 P 

Then the variation in entropy becomes 

d 4 p 



(p°)6(p 2 +M 2 )p»p 6f(p) = 



5S = -2tt / d A X 



(27T) 



I 0(p°),5(p 2 +M 2 )p 



[1 + fea (p)] /e 9 (p) 



(<V) 2 



(11.17) 



(11.18) 



In the classical theory, the distribution function is concentrated on the positive frequency mass shell. 
Therefore, it is convenient to label momenta just by its spatial components p, the temporal component being 
necessarily ui p — y/ M 2 + p 1 > 0. In the same way, it is simplest to regard the distribution function as a 
function of the three momentum p alone, according to the rule 



/< s > (x,p) = /[x, K,p)] 



(11.19) 



where / represents the distribution function as a function on four dimensional momentum space, and /*- 3 ' its 
restriction to three dimensional mass shell. With this understood, we shall henceforth drop the superscript, 
using the same symbol / for both functions, since only the distribution function on mass shell enters into 
our discussion. The variation of the entropy now reads 



ss = -l fd'x f d3p 1 

2i i (2nf{l + fe (p)]fe q (p) 



(SfY 



(11.20) 



From Einstein's formula, we conclude that, in equilibrium, the distribution function is subject to Gaussian 
fluctuations, with equal time mean square value 



Sf (t, X,p) Sf (t, Y, «f) ) = (2tt) 3 S(X-?)S (p- q) [1 + f eq (p)] f eq (p) 



(11.21) 



One of the goals of this paper is to rederive this result as the kinetic theory limit of the general fluctuation 
formula given for the propagators in the Introduction, Eq. (1.6). For the time being, we only observe that 
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this fluctuation formula is quite independent of the processes which sustain equilibrium; in particular, it 
holds equally for a free and an interacting gas, since it contains no coupling constants. 

In the interacting case, however, a stochastic source is necessary to sustain these fluctuations. Following 
the discussion of the FDR above, we compute these sources by writing the dissipative part of the equations 
of motion in terms of the thermodynamic forces 



F(X,p) 



1 



Sf (X,p) 



[l + fe q ( P )]fe q (p) (27T) 



(11.22) 



To obtain an equation of motion for / (X,p) multiply both sides of the Boltzmann equation Eq. ( [I.lOj ) by 
6 (fc°) S (k 2 + M 2 ) and integrate over k° to get 



df 
dt 



-v/ 



Lcol 



Upon variation we get 



d(Sf) k 1 
— ^— + — V(d/J = — dlcoi 
at out <^k 



(11.23) 



(11.24) 



When we write 5I co l in terms of the thermodynamic forces, we find local terms proportional to F (k) as 
well as nonlocal terms where F is evaluated elsewhere. We shall keep only the former, as it is usually done 
in deriving the "collision time approximation" to the Boltzmann equation E(| (also related to the Krook - 
Bhatnager - Gross kinetic equation), thus we write 



SIcoi (k) ~ -u k is 2 {X, k)F(X, k) 



where 

u 2 (X, k) 
fc° = u> k , and 



d A Pi 



Ha 



(p°)5{p 2 t -M 2 



(27T) 5 (pi + p 2 - P3 - k) 



[1 + f eq (px)] [1 + f eq (pa)] f eq (ps) f eq (k) 



(11.25) 



(11.26) 



(11.27) 



Among other things, the linearized form of the Boltzmann equation provides a quick estimate of the 
relevant relax ation time. Let us assume the high temperature limit, where / ~ T/M, and the integrals 
in Eq. QEI.26P are restricted to the range p < M. Then simple dimensional analysis yields the estimate 
t ~ M I X 2 T 2 for the relaxation time appropriate to long wavelength modes . 



C. Fluctuations in the Boltzmann equation 



Obs ervance of the FDT demands that a stochastic source j be present in the Boltzmann equation Eq. 
( [I.10| ) (and its linearized form, Eq. ( 1I.24|) ) which should assume the Langevin form: 



% + —Vf=—I col+J (X,k) 
at uj k uj k 



Then 



(j(X,p)j(Y,q))=- 



1 6I col (X,p) 1 5I co i{Y,q) 



lo p 5F{Y,q) Lu q 5F{X,p) 



(11.28) 



(11.29) 
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From Eqs. @L2|), |L2|) and |L27]) we find the noise auto-correlation 



j(X, k)j {Y,p)) = 26^ (X-Y)s(k- p) v 2 {X, k) 



(11.30) 



where v 2 is given in Eq. ( !L2q ). Eqs. (11.30) and (11.26) are the solution to our problem, that is, they 
describe th e fluct uations in the Boltzmann equation, required by consistency with the FDT. Observe that, 
unlike Eq. (11.21), the mean square value of the stochastic force vanishes for a free gas. 

In this discussion, of course, we accepted the Boltzmann equation as given without tracing its origin. 
We now want to see how the noises in Eq. ([1.30) originate from a deeper level, that related to the higher 
correlation functions, which we call the correlation noises. 



III. KINETIC FIELD THEORY, FROM DYSON TO BOLTZMANN 

Our goal in this section is to show how the Boltzmann equation arises as a description of the dynamics of 
quasiparticles in the kinetic limit of field theory. To this end, we shall adopt the view that the main element 
in the description of a noncquilibrium quantum field is its Green functions, whose dynamics is given by the 
Dyson equations. This connects with the results of our earlier paper on dissipation in Boltzmann equations 
||10|| . The task is to find the noise or fluctuation terms. The need to upgrade the Boltzmann equation to a 
Langevin form will lead to a similar generalization of Dyson's equations, whose physical origin will be the 
subject of the remaining of the paper. 

The discussion of propagators is simplest for a free field theory, and so, following our choice of physical 
clarity over formal rigor in the exposition, we shall first discuss nonequilibrium free fields. The general case 
follows. 



A. Free fields and propagators 

Let us focus on the nonequilibrium dynamics of a real scalar quantum (Heisenberg) field $ ( x) , obeying 
the Klein - Gordon equation 

(□ - m 2 ) $ (x) = (III.l) 
and the canonical equal time commutation relations 

$ {x , t) , $ (y, *)] = -ihS (x - y) {111.2) 

(from here on, we take U = 1). 

We shall assume throughout that the expectation value of the field vanishes. Thus the simplest nontrivial 
description of the dynamics will be in terms of the two - point or Green functions, namely the expectation 
values of various products of two field operators. Of particular relevance is the Jordan propagator 

G(x,x') = ([^(x)^(x')]) (ni-3) 

which for a free field is independent of the state of the field. From the Jordan propagator we derive the 
causal propagators, advanced and retarded 

G adv (x,x') = -iG(x,x') 9 (t'-t), G ret (x,x') = iG(x,x')9(t-t') (IH.4) 

These propagators describe the evolution of small perturbations (they are fundamental solutions to the 
Klein - Gordon equation) but contain no information about the state. For that purpose we require other 
propagators, such as the positive and negative frequency ones 

G+ (x, x 1 ) = ($ (x) $ (a;')) , G_(x,a/) = ($(x')$(x)) (III.5) 
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Observe that G = G+ — G- . The symmetric combination gives the Hadamard propagator 

G 1 =G++G-=({$(x),$(x')}) (IH.6) 

Note that while the Jordan, advanced and retarded propagators emphasize the dynamics, and the negative, 
positive frequency and Hadamard propagators emphasize the statistical aspects, two other propagators 
contain both kinds of information. They are the Feynman and Dyson propagators 

G F (x, x') = (T [$ (x) $ (x')}) = i [Gi (x, x') + G (x, x') sign (t ~ t')] (111.7) 

G D (x, x') = (f [$ (x) $ (x')} ) = g [ G i ( x ' X ')~ G ( x > x ') si 9 n (* " 0] ( IIL8 ) 
where T stands for time-ordered product 

T [$ (x) $ (a/)] = $ (x) $ (x') d(t-t') + ® (x') $ (x) 9 (f' - t) (HI.9) 
and T for anti- temporal ordering 

f [$ (x) $ (a/)] = $ (a;') $ (a;) (i - *') + $ 0*0 $ 0') (*' - t) (IIL 10) 



B. Equilibrium structure of propagators 



In this subsection, we shall re view s everal important properties of the equilibrium propagators which follow 
from the KMS condition [Eq. ( III. 1 2 1 below) f47|j , and general invariance properties. 

In equilibrium, all propagators must be time-translation invariant, and may be Fourier transformed 



G(x,x') 



d 4 k 



e ik ( x - x ')G(k) 



(III.ll) 



In particular, because the Jordan propagator is antisymmetric, we must have G (cj, kj = — G u, kj . Also, 

since G(x,x') =G(x',x)* = —G(x,x')*, G (k) =G(k)* . 

The positive and negative frequency propagators are further related by the KMS condition 



G+ [(t, x) , (f , x')] = G_ [(t + 0, x) , (f , x 1 )] 
where j3 is the inverse temperature. With G+ — G_ = G, we get 



G+ (k) 



G(k) 
1 - e^ 



sign (k 



(*°) 



1 

s /3|fe°| _ 1 



G(fc) 



(III. 12) 



(III. 13) 



Adding these two equations, we find 

Gi (k) = 2sign (k°) 



1 



e P\k°\ - 1 



G(fc) 



1 



s /3|fc°l _ 1 



G(fc) 



(III. 14) 



(III. 15) 



We may co nsider this formu la as the quantum generalization of the FDT, as we shall see below. Let us stress 
that Eqs. (111.12) to (111.15) hold for interacting as well as free fields. 
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Of course, for the homogeneous solutions to the Klein - Gordon equation (G, G + , G_ and G\) we must 
have 



G{k) = 6(k 2 + m 2 ) g (k) 



(111.16) 



which leads to 



G re t (x, x') 



d 4 k 



(27rr-(fco-fer+^ 



g[uk,k 



2tt 



(111.17) 



and to 



Gf {x, x') 



d^k 



H) 



2ttS (k 02 - u>l) 



-k 02 + u>i — ie 



e /3|fe°l _ i 



2tt 



(111.18) 



with similar formulae for G ac iv and Gd, respectively. It is remarkable that all propagators may be split into 
a vacuum and a thermal contribution, with the thermal part being the same for all propagators except G, 
G re t and G a d v , where it vanishes. Also, we have expressed all propagators in terms of g; in the language of 
the Lehmann decomposition, this is just the density of states p8[ . 

We shall finish this subsection by expanding our remark on Eq. (111.15) being the fluctuation dissipation 
theorem p9[ . Suppose we try to explain the quantum and statistical fluctuations of the field by adding an 
external source -j (x) to the right hand side of the Klein - Gordon equation ( III.l ). The resulting field would 
be 



$ (x) — \ d X Gret (x, x) j {x) 



If the process is stationary 



we get 



(J(x)j(x')) 



d 4 k 



(111.19) 



u{k) 



Gi (k) 



2\G ret (fc)r 



From Eqs. (lTLT5|) and (111.17) 



y(k) 



e P\k°\ - 1 



\lmG-e\{k)\ 



(111.20) 



which is a generalized form of the FDT, including both quantum and thermal fluctuations. 

So far, we have intentionally left everything expressed in terms of the density of states g (k). For a free 
field, we can compute this explicitly 



g (k) = 2nsign (fc°) 
with which we can fill in the remaining results. 



(111.21) 
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C. Interacting fields and the Dyson equation 



Let us now consider a weakly interacting field, obeying the Heisenberg equation 



□ 



$ (x) - (a) = 



(111.22) 



and the same equal time canonical commutation relations, Eq. (III. 2). As before, we shall assume that 
the expectation value of the field vanishes identically, and seek to describe the dynamics in terms of the 
propagators introduced earlier. 

In the usual approach to field theory, where one focuses on computing the S-matrix elements, rather than 
the causal evolution of fields, the leading role is played by the Feynman propagator, which is directly related 
to the S-matrix through the LSZ reduction formulae, and has a simple perturbative expansi on E8| , |50| , |5T[ . 
We may obtain a dynamical equation for the Feynman propagator by noting that, from Eq. ( [II. 9| ) 



(-□ + m 2 ) T [$ (x) $ (x')] = T [(-□ + m 2 ) $ (x) $ (x')] - i5 (t - t') 



Therefore 



(-□ + m 2 ) G F (x, x') = -id (t - t') - ~ (T [$ 3 [x) $ (:r/)] ) 



(111.23) 



(cfr. Eq. (III. 18)). This is the Dyson equation for the propagator, relating the evolution of the Feynman 
propagator to higher order (in this case, four point) correlation functions. 

As different from an IN-OUT matrix element of the S-matrix, in this case we have an IN - IN expectation 
value taken with respect to a nontrivial state defined at some initial time. Thus the perturbative expansion 
of the self energy term cannot be expressed in terms of the IN - IN Feynman propagator alone. We should 
rather have 



(T [$ 3 (x) $ (a:')] ) ~ 3G F (x, x) G F (x, x') - i\ / d 4 y {G 3 F (x, y) G F (y, x') ~ G 3 _ (x, y) G+ (y, x')} 



(since we use full propagators in the internal lines, two-particle reducible (2PR) graphs must not be included). 
Thus to obtain a self-consistent dynamics, we must enlarge the set to include other propagators as well. Of 
course, we are assuming that the initial state is such that Wick's theorem holds (for example, that it is 
Gaussian) - this issue is discussed in detail in [^0) . We are also leaving aside issues of renormalization [[52| . 

We overcome this difficulty by adopting as fundamental object the closed- time-path ordered propagator 
Gp (x a ,y b ). This object is equivalent to four ordinary propagators: if we write Gp (x a ,y b ) — G ab (x,y), 
then G 11 (x, y) = G F (x, y), G 12 (x, y) = G_ (x, y), G 21 (x, y) = G+ (x, y) and G 22 (x, y) = G D (x, y) . 

We can obtain closed dynamical equations for these four propagators. Actually there is a slight redundancy, 
but this set has the advantage of being very simple to handle (see Chou et al in [^9[ for details) . The equations 
read 



-□ + m 2 + ^-G F (x,x) 



G ab (x,x') 



i\ 2 



c cd / d 4 y S ac (x, y)G db (y, x') = -ic ab 5 (x - x 1 ) 



(111.24) 



Z ac (x,y) = [G ac (x,y)} 3 (111.25) 

The matrix c(cn = c 11 = 1, C22 = c 22 = —1, all others zero) keeps track of the sign inversions associated 
with the reverse temporal ordering of the second branch. This form of the Dyson equation is relevant to our 
discussion. 
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D. The kinetic theory limit 



In equilibrium the propagators are time-translation invariant. Out of equilibrium this is no longer true. In 
the kinetic theory regime, however, the propagators depend mostly on the difference variable u = x — x', with 
the corresponding Fourier transform depending weakly on the center of mass variable A — (1/2) (a; + x'). 
As such, the propagators take the form 



74 7 

G ab {x,x') = I ie lk ( x - x ')G ab (X,k) 



(27T) 



The E kernel has a similar expression 



T, ab (x,x') = I -^^e lk ( x - x ')^ ab (X,k) 



(111.26) 



£° b (A,fc) = J jQ l^ G « b {X, Pi )\ [(!fcr) 4 f (J> * 



The weak dependence on X allows for the approximations (details in |icj) 

G ab (x,x) = [ ^G ab (x,k) r 1 ' !>k 
J (2tt) 4 



(27T) 



jG ab {X,k) 



(111.27) 



d*y Z ac (x,y)G dh (y,x') 



_d*k_ 
(2^ 



ik(x—x'\ 



E ac (X, k) G db (X, k) 



and the equations of motion become 



k 2 - ik^ 



d 1 



dX^ 4 



-a x + M' (X) 



G ab (A, k) - —c cd £ ac (A, k)G db (A, k) = -ic ab 



(111.28) 



M 2 (X) = m 2 + X f ^ <-" lh 



2 7 (27r)' 



G a0 {X, k) 



(111.29) 



Alternatively, we may think of the propagators as functions of x' , leading to an equation of the form (cfr. 
Eq. (pL24|)) 



-□' + m 2 + ^G F (x',x') 



In the kinetic limit, this yields 



G ab (x, x') - l -^c cd J d 4 y G ac (x, y)H db (y, x') = -ic ab 5 (x - x') (111.30) 



k 2 +ik fi - 



d 1 



-U X +M 2 (X) 



G ab (X, k) - ^-c C(i G ac (X, k)Y, db {X, k) 



dX» 4 

Taking the average and the difference of Eqs. ( 111.28 ) and (|lll.3l| ) we get 



-ic 



ah 



(111.31) 



k 2 ~-a x + M 2 (A) 



7 A 2 

G ab {X, k) - —c cd {E ac (A, k)G db (A, k) + G ac (X, k)H db (A, k)} = -ic ab (111.32) 



kvJL.G ab {X, k) + ^c cd {S QC (A, k)G db (A, fc) - G ac (A, k)Y, db (A, k)} = 



(111.33) 
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We recognize the first equation as a mass shell condition on the noncquilibrium propagator. The second 
equation is the kinetic equation proper, describing relaxation towards equilibrium. 

To investigate further this equation, we observe that since both terms are already of second order in A (see 
Jul), it is enough to solve the mass shell condition to zeroth order. That is, we assume that the renormalized 
mass M 2 is actually position independent, and write 



G ab (X, k) = G% b (M 2 ,k) + Gf ta 



(X,k) 



(111.34) 



where the Gq (M , k) are vacuum propagators for a free field with mass M , and G s t a t is the non vacuum 
part 



Gf tat (X,k) = 2TrS{k 2 



+ M 2 ) f (X, k) 



(111.35) 



which we assume is the same for all propagators involved, as in the free field case. / (X, k) has the physical 
interpretation of a one p articl e distr ibutio n func tion for quasi particles built out of the field excitations. 
Substituting Eqs. ([11.34) and (III. 35) into (111.33), and assuming, for example, that k° > (/ must be even 
in k, because of the symmetries of the propagators) immediately shows that the dynamics of / is given by 
the Boltzmann equation flll.lOD , ( |EI. 1 1[ ) and ( Jll. 12| ) . 

We shall not discuss further the region of validity of the hypothesis underlying the kinetic limit, except 
to observe that this issue is far from trivial. On general grounds, one expects that propagators will depend 
strongly on the difference variable on scales tc ~ M -1 . For smooth initial conditions, the scale for depen- 
dence on the average variable is set by the relaxation time r ~ M/X 2 T 2 . A nontrivial kinetic limit exists if 
t 3> tc- Already this simple estimate shows that one would expect trouble in theories with strictly massless 
particles, such as gauge or Goldstone bosons jjl2|. If particle masses are not specially protected, then at large 
temperature the physical mass M ~ vAT, and tc/t ~ A will in general be suitably small. 



E. Stochastic Dyson equations 

Our derivation of kinetic theory from the Dyson equations leads to an incomplete Boltzmann equation: 
in addition to the usual collision integral an explicitly stochastic term ought to be in place. Mor eover , 
fluct uation - dissipation considerations demand that this stochastic term has auto-correlation Eqs. ( II. 3C ) 
and (II. 26). Since no manipulation of the deterministic Dyson equations will yield a stochastic term like 
this, we posit that when quantum field theory is viewed in the statistical mechanical context, the Dyson 
equations themselves are incomplete. Suppose we add a stochastic driving term F ab to them (we shall justify 
this later): 



-□ + m 2 + -G F (x,x) 



-W + m 2 + ^G F (x',x') 



G ab (x, x') - l -^-c cd J d 4 y S QC (x, y)G db {y, x') = -ic ab 6 (x - x') - iF ab (x, x') 



(111.36) 



G ab (x, x') - ^c cd I d 4 y G ac (x, y)Y, db (y, x') = -ic ab 5 (x - x 1 ) - iF ab {x, x') 



(111.37) 



In the kinetic limit, the random forces become 



F ab (x,x r ) = 



d*k 

(£7 



e tk(x-x')pab fc ) 



(111.38) 



(and similarly for F). Leaving aside the random fluctuations of the mass shell, we find the new kinetic 
equation 
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fc"- 



-G ab (X, k) + —c cd {E QC (Af, k)G dh (X, k) - G ac (X, k)Y, db (X, k)\ = H ab (X, k) 



(111.39) 



where 



l r 

2 



F — F 



(X,k) 



(111.40) 



Our problem now is to justify changing the Dyson equation to Eq. (|IL36|), and to expound the physical 
meaning of this new stochastic equation. To do this we need to use functional methods, to which we now 
turn. 



IV. CORRELATION NOISE AND STOCHASTIC BOLTZMANN EQUATION 

Our goal in this section is to show how noise terms such as those introduced above from phcnomenological 
considerations may actually be systematically identified from an appropriate effective action. In this section 
we shall limit ourselves to finding a suitable recipe to identify the noise terms, and to compare the results 
to the phenomenological discussion above. The physical foundations of the recipe shall be discussed in the 
following section. 



A. Fluctuations in the propagators 



We shall now adapt the foregoing discussion to the study of fluctuations in the dynamics of the two point 
functions. The first step is to notice that this dynamics can be obtained from the variation of the 2PI action 
functional (we derive this formula in Section V, Eq. ( V.52)) 



1 



r [G ab ] = — In [DetG] - -c ab / d 4 x (-□ + to 2 ) G ab {x, x) 



(IV.l) 



A 



Cabcd / d*xG ab (x,x)G ab (x,x) 



+ l -^c abod c efgh I d 4 xd 4 x' G ae (x, x') G bf (x, x') G cs (x, x') G dh (x, x') 



The resulting equations of motion 
1 



—r- 1 - 
2 ab 2 



Cab (-□ + m 2 ) + -c abc dG cd (x, x) 



i\ 2 3 
6 (x, x') + —c ac Cbd [G cd (x, x')] ' = 



(IV.2) 



are seen to be equivalent to the Dyson equations Eq. ( |III.24| ). 

As we discussed in the Introduction, we shall incorporate quantum fluctuations in the evolution of th e 
Green function G ab by explicitly adding a stochastic source {—l/2)n a b to the right hand side of Eq. ( [V.2| ). 
Let us write G afc = G a b + A ab , and expand the 2PI CTP EA to second order (the first order term vanishes 
by virtue of Eq. fliV.2|)), 



ST = s 2 r 



(IV.3) 



6 2 T [A ab ] = ^G~ b 1 A bc G~ d 1 A da - ±c abcd J d 4 x A ab (x, x) A ab (x, x) 



(IV.4) 



i\ 2 



+^-c abcd c efgh / d A xd A x' G ae (x, x') G bf (x, x') A c ^ (x, x') A dh (x, x') 
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From now on, wc shall assume that the background tadpole vanishes, and identify the mass with its renor- 
malized value. We now have, as discussed in the Introduction 



<A ah A crf ) = i 



S 2 T 



_6G ab SG cd _ 

To sustain these fluctuations, the noise auto-correlation must be 



That is 



(KabKcd) = (4«) 



SA ab SA cd 



5 2 r 



1 1 



A=0. 



(n ab (x, x') K cd (y, y')) = N abcd (x, x' , y, y') + N l a f cd (x, x',y, y') 



(IV.5) 



(IV.6) 



Nabcd (x, x', y, y') = G d l (y\ x) G 6 / (x , y) + Gj (x, y) G db (y\ x) 

Nltd (x, x', y, y') = X 2 [G eg G fh ] (x, x') c acef c bdgh S (x-y)S (x' - y') 
-iXcabcdS (x - x') 5(y - y') 



B. Free fields 



Let us begin by askin g whe ther for free fields the quantum fluctuations Eq. (1V.5) go into anything like 
the classical result Eq. (11.21) in the kinetic theory limit. There is no obvious reason why it should be so, 
s ince the physical basis for either formula is at first sight totally different. As we saw in the Introduction, Eq. 
(1V.5) simply reproduces the full quantum fluctuations, computed in terms of the propagators themselves on 
the assumption that Wick theorem holds (which is an assumption on the allowed initial states of the field, 
sec p|) 



(A ab A cd ) 



5 2 T 



SG ab SG cd 



QacQdb _|_ £ida£ibc 



(IV.7) 



while the classical auto-corr elatio n Eq. ( 11.21 ) has been found by applying Einstein's formula to the phe- 
nomcnological entropy eq. (11.13). The only clear point of contact between both approaches is that both 
assume Bose statistics. 

Introducing the Wigner transform of the fluctuations 



A ab (x,x') 



tk( X - X ') A ab k y 



X = \{x + x'). 



(IV. 



we observe that in this case we are not entitled to assume that the dependence of the Wigner transform on 
X is weak. Eq. (IV.8|) has a formal inverse 



A ab {X, k) = J du e iku A ab (x + |, X 



(IV.9) 



and from Eq. (IV.7) we get 

<A ab (X,p) A cd (Y, q)) = I dudv e i( - pu+qv >K [V, Y, u, v] , 
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where 



K[X,Y,u 1 v] = G ac [X+^Y + -)G bd [X-^,Y--) + G ad [X + ^,Y - -) G bc [X - -,Y , 
The propagators in the right hand side are equilibrium ones, and so we can use the representation Eq. 

(mil) 



K[X,Y,u,v] = / ^- —t / dudv e l{pu+qv) K\X,Y,r,s} 

J (2tt) 4 (2tt) 4 J 



K[X,Y,r,s] 



Now integrate over u, v and s 



(A ab (X,p)A cd {Y,q)) 



e ir(X -Y+±(u-v)) e is(X -Y - ±(u-v)) Qac / \ gr&d / \ 
+e tr(X-Y+±(u+v)) e ts(X-Y-±(u+v)) G ad / s g&c / s 



16 / d 4 r e '2(r+p)(x-y) (IV1Q) 
[<5 (p + g) G ac (r) G M (r + 2p) + 6 (p - g) G ad (r) G bc (r + 2p)] 



We have 16 different quantum auto-correlations to compare against a single classical result, so we can only 
expect real agreement in the large occupation number limit, where all propagators converge to the same 
expression. With this proviso in mind, we can choose any combination of indices to continue the calculation. 
The most straightforward choice (to a certain extent suggested by the structure of the closed time - path, see 
p6| , p7| ) is a = b = 1, c = d = 2; we are thus seeking the correlations among the fluctuations in the Feynman 
and Dyson propagators 

(A 11 (X,p) A 22 (Y, q)) = 16 (2ir) 2 [S (p + q) + S (p - g)] J d 4 r e ^(r+p)(x-Y) s p + M 2] g ^ + 2p) 2 + Af2 

[0 (-r°) + f eq (r)] (-r° - 2 P °) + / eg (r + 2p)] 
The arguments of the delta functions can be simplified 

(IV.11) 



(A 11 (X,p) A 22 {Y, q)) = 4 (2tt) 2 [5(p + q) + 5(p- q)} I d 4 r e ^+p)(x-Y) 6 p + M 2j 

J [rp + p 2 ] [0 (-r°) + f eq (r)] [0 (-r° - 2p°) + / e , (r + 2p)] 

A difference from the classical case already stands out here: in the quantum case, a fluctuation in the 
number of particles with momentum p correlates not only with itself, but also with the corresponding 
fluctuation in the number of antiparticles with momentum —p. This is unavoidable, given the symmetries 
of the propagators in this theory. 

Let us stress that we are trying to push the quasi - particle (kinetic) description of quantum field dynamics 
beyond the calculation of mean values (of such quantities as particle number or energy density), to account 
for their fluctuations. The calculation of the fluctuations of the distribution function for on-shell particles 
gives a crucial consistency check on suc h an a tt empt. Indee d, we k now that each on-shell mode of the free 
field contributes an amount [cfr. Eqs. ([11.15), (IV. 12) and ( III.2l | )1 pk ~ <jJk (1/2 + / e <j)to the mean energy 
density, where f eq is the equilibrium distribution function Eq~ (f[I.15| ). The fluctuations of this quantity 
at equilibrium will be given by (Sp 2 ) — T 2 (dpk/dT) ~ u 2 f eq . (1 + f eq ). So, if these fluctuations are still 
described by a distribution function consistent with ordinary statistical mechanics, then this distribution 
function must fluctuate like in Eq. (11.21). (This may in the face of it be a rather big if). 

For large M 2 , the condition that p is nearly on - shell means that the spatial components are much smaller 
than the time component, and we may approximate 
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S [r 2 + M 2 ] 5 [rp + p 2 ] ~ 6 [p 2 + M 2 } —8 (r° + p°) 



thus obtaining 

(A 11 (X,p) A 22 (Y, q)) | on _ shoU = ±- (2nf [S( P + q ) + S(p- q)\ f eq (p) (1 + f eq (p)) J [p 2 + M 2 ] s(x-y) 
To finish the comparison, assume, e. g., that p° > 0, then 



S(p-q)S [p 2 + M 2 ] =S(q°- uo q ) S(q-p)S [p 2 + M 2 ] 

= 2cj q 9 (q°) 6 (q- p)S[p 2 + M 2 ] 6 [q 2 + M 2 ] 

This results suggests writing 

A 11 (X,p) = 2nSf {X,p)5 [p 2 + M 2 ] + off - shell terms. 
Taking p° and q° to be positive, this yields 

[Sf (t, X,p) Sf (t, Y, q) ) = (2tt) 3 S (q- p)5(X-?) f eq (p) (1 + f eq (p)) 



(IV.12) 



(IV. 13) 



(IV. 14) 



which is identical to Eq. (11.21). This is one of the most i mportant results of this paper, as it gives a whole 
new meaning to the phenomenological entropy Eq. (11.13) 

We have thus completed our proof, and obtained new independent confirmation of the validity of our 
scheme for introducing fluctuations in the dynamics of correlations. 



C. Interacting fields and the Boltzmann - Langevin equation 

The results of the previous section already imply that the full stochastic Dyson equation will go over to 
the Boltzmann - Langevin equation in the kinetic limit. Indeed, the structure of the fluctuations does not 
change drastically when interactions are switch on, and since they become identical in the classical limit, 
the noise in the Dyson equation necessary to sustain the fluctuations at the quantum level must go over to 
the noise in the Boltzmann equation, which plays the same role in the classical theory. Nevertheless, it is 
worth identifying exactly which part of the quantum source auto - correlation goes into the classical one in 
the correspondence limit. 

Concretely, our aim is to begin with the stochastic Schwinger - Dyson equation 



-G 



a b 



Cab (-□ + m 2 ) + ^c a6cd G cd (x, x) 



iX 2 

S (x, x 1 ) + —c ac c bd [G cd (x, x' 



1 



n ab (IV.15) 



(111.36) and (111.37) 



where the noise auto - correlation is given by Eq. (IV.6). We then identify the forces appearing in Eqs 



F ab (x, x')=i d A y c ac n cd (x, y) G db (y, x') 



(IV.16) 



f\ „db 



F ab (x, x') = i d 4 y G ac (x, y) K cd (y, x 1 ) ti 



In condensed notation, 



F -F 



f ac fidb s~iac db\ 
{c K cd G - G K cd C ) 



(IV. 17) 
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whose Wigner transform plays the role of random force in the kinetic equation (111.39). Restricting ourselves 
to on-shell fluctuations, we can compute the auto-correlation of this force and compare the result to the 
classical expectation Eq. ( 1 1 . 3 C| ) . 

Let us observe from the outset that the classical result involves the expression v (Eq. flH.26|)), w hich, 
through I + (Eq.( |II.27 )), is related to the Fourier transform of the cube of a propagator. In Eq. (IV.6), the 
first term N contains the inverse propagators, which in turn is related to the cube of propagators through 
the Dyson equations ( IV.2j ). The other term N lnt contains no such thing. Thus it is clear that our only 
chance lies in the first term, the other one contributing to sustaining the nonclassical correlations already 
present in the free field case. Correspondingly, we shall ignore N lnt in what follows. 

We thus approximate 



(KabKcd) — G. G 



be 



G-l /~f— 1 
ac U db 



(IV. 18) 



leading to 



(Jjabjjcdj _ [ c ae G /b _ G ae c fb^ ^cg Q hd „ Qcg ^ G]£Gj* + G^G^j 



= - \-c ae G fb + G ae c fb ] 
4 L J 



cd/^t—lc , c-d/^r— lc £C s~r— Id cc/^—ld 
°e G f + d f G e ~0 f G e ~0 e G f 



= - [G ad G 



ad/~t—lbc 



Q—lad^tbc QbdQ—lac Q—lbd qcic _j_ ^ ^ ac ^bd c^C^^J 



For the same reasons as above, we shal l disregard the propagator - independent terms. 
Next, we write (recalling Eq. ( III. 25 )) 



H ab ( Xl x') 



(rv.19) 



H ab (X, k) = J d 4 u e~ lku H ab (X + H X 



to get 



(H ab (X,p) H cd (Y, <l)) = \j d4 ud 4 v e - l(pu+qv) [J - K] 



(IV.20) 



where, using the translation invariance of the equilibrium propagators 



J = 



d A r d s 
^77T ~J exp 



(27T)* (27T) 



ir X -Y 



u + v 



{G ad (r) G- lbc (s) + G- lad (r) G bc (s)} 



exp 



is \X-Y- 



u + v 



K = 



d A r d Is 
j j ex P 

(27T) 4 (27T) 4 



ir X -Y - 



{G bd (r) G- lac (s) + G- Ud (r) G ac (s)} 



exp 



is X -Y 



Upon integration over u and v, the K term gives a contribution proportional to 5^ (p + q). This is 
unrelated to the noise auto - correlation, being only a cross correlation between the positive and negative 
frequency components of the source, and we shall not analyze it further. We also restrict ourselves to the 
case where a = b = 1, c = d = 2. Using the reflection symmetry of the equilibrium propagators, obtaining 
the inverse propagators from Eq. ( IV. 2 ) and retaining only the dominant term in the correspondence limit, 
we find 



19 



4A 2 f 

(H 11 (X,p)H 22 (Y,q))~—5(p-q) J d A s cos[2{ S +p)(X-Y)]V 12 {s + 2p)G 12 { S ). (IV.21) 

An analysis of this expression shows that in the high temperature limit the correlation length is of order 
M _1 . This is a microscopic scale, much smaller than the macro scales of relevance to the kinetic limit (if 
this limit exists). Therefore we are justified in writing 

(H 11 (X,p) H 22 (Y, q)) ~ 7,5 (X - Y) (IV.22) 



We compute 7 by simply integrating Eq. (IV.21) over X 



7 = - (2ny £ (p) G (p) 5(p-q). (IV.23) 



From Eqs. (|III.27D and ( |[I.26|) we get 

7 = w p (2tt) 2 6{p-q)6 (p 2 + M 2 ) v 2 (IV.24) 

Assuming p°, q° > 0, this is 

7 = 2w 2 (2tt) 2 5 (p - q) S (q 2 + M 2 ) 5 (p 2 + M 2 ) v 2 (IV.25) 

So, writing 

H (X, k) = 2iruj k 6 (k 2 + M 2 ) j (x, k\ + off - shell (IV.26) 

we get the final result 

(j (X,p)j (Y, q)) = 28 (p- q)6(X~ Y) v 2 , (IV. .27) 



which agrees with the classical result, Eq. ( II.30| ) 



We have shown that there is a piece of the full quantum noise which can be identified with the classical 
source j. Clearly j does not account for the full quantum noise, the difference being due among other things 
to the role of negative frequency in the quantum theory. 

Finally, we note that Abe et al. |l5| have given a nonrelativistic derivation of the Boltzmann - Langevin 
equation, while ours is fully relativistic, being also immune to the reservations expressed by Greiner and 
Leupold (T|. 

V. MASTER EFFECTIVE ACTION 

So far in the paper we have referred several times to the possibility of conceiving the low order correlations 
of a quantum field as the field variables of an open quantum system, interacting with the environment 
provided by the higher correlations. The goal of this final section is to present a formalism, the master 
effective action, built on this perspective. In particular, in this formalism the usual Dyson equations are seen 
to emerge from the averaging over higher correlations. As a simple consequence and illustration, we derive 



Eq. (IV.l) 



A. The low order effective actions 

The simplest application of functional methods in quantum field theory concerns the dynamics of the 
expectation value of the field f33| . The expectation value or mean field may be deduced from the generating 
functional W [J] 
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exp{iW[J}} = J £>$ expjiS'^+t J d 4 x J (x) $ (x)\ (V.l) 

° J J=0 

We obtain the dynamics from the effective action, which is the Legendre transform of W 

r [0] = W [J] - J d 4 xJ (x) cj) (x) (V.3) 
The physical equation of motion is 

In a causal theory, we must adopt Schwinger's CTP formalism. The point x may therefore lie on either 
branch of the closed time path, or equivalently we may have two background fields cj) a (x) — <fi(x a ). The 
classical action is defined as 

5[$ a ] = Sft 1 ] -S[<£ 2 ]* (V.5) 
which automatically accounts for all sign reversals. We also have two sources 

and obtain two equations of motion 

W = ° <™> 

However, these equations always admit a solution where (f) 1 = 4> 2 — 4> is the physical mean field, and after 
this identification, they become a real and causal equation of motion for <p. 

The functional methods we have used so far to derive the dynamics of the mean field may be adapted 
to investigate more general operators. In order to find the equations of motion for two-point functions, for 
example, we add a nonlocal source K ab (x, x') pl]JlO|| 



exp{iW[J a ,K ab }} = J expi js[$ a ]+ J d i xJ a <i> a + ^ J d i xd A x' K ab <P a <P b 



(V.7) 



It follows that 



-= l -[^{x)<p b {x') + G ab (x,x') 



5K ab (x, x') 2 

Therefore the Legendre transform, the so-called 2PI effective action 



T [r, G ab ] = W [J a , K ab ] - J d A x J a ct> a - ^ J d^xd^x 1 K ab [r<f> b + G ab ] (V.8) 
generates the equations of motion 

^r = -^-W;^ = ~^ (v.9) 

The goal of this section is to show these two examples as just successive truncations of a single object, the 
master effective action. 
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B. Formal construction 



In this section, we shall proceed with the formal construction of the master effective action, a functional 
of the whole string of Green functions of a field theory whose variation generates the Dyson - Schwinger 
hierarchy. Since we are using Schwinger - Keldish techniques, all fields are to be defined on a closed time 
path. Also we adopt DeWitt's condensed notation [ p4j . 

We consider then a scalar field theory whose action 

S[$] = \s 2 & + S int [$] (V.10) 
decomposes into a free part and an interaction part 

oo 1 

s mt m = Y,ri s ^ n (v - n) 

71 — 3 

Here and after, we use the shorthand 

K n <S> n = J d d x 1 ..J d x n K nal ... an (x 1 ,...x n )$ a \x 1 )...$ an (x n ) (V.12) 

where the kernel K is assumed to be totally symmetric. 
Let us define also the 'source action' 

J[$] = Jl$+i/ 2 $ 2 + JinM (V.13) 

where Ji n t [$] contains the higher order sources 

JinM = Y] ~. J»# B (V.14) 

and define the generating functional 

Z[{J n }\ = e'WK-M] = / D$ e *St[»,{J»}] (V.15) 



m 

n— 3 



where 

St[$, {J n }} = Ji$ + i(5 2 + J 2 )$ 2 + S int [$] + J int [$] (V.16) 

We shall also call 

S int [$] + J int [$]=S I (V.17) 

As it is well known, the Taylor expansion of Z with respect to J\ generates the expectation values of path 
- ordered products of fields 

tj 7 r Z XT 7 v = (P{* a \x 1 )...** n (x n )})=F?"-* n (xi,...x n ) (V.18) 
6J la i{xi)...bJi an [x n ) 

while the Taylor expansion of W generates the 'connected' Green functions ('linked cluster theorem' Q) 
S n W 



SJ la i(xi)...SJ 1 



= (P{$° ( Xl )...$ a (x„)}) conncctcd = CI - a (xi, ...x„) (V.19) 
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Comparing these last two equations, we find the rule connecting the f s with the C's. First, we must 
decompose the ordered index set (ii,..A n ) (ik — {xk,a k )) into all possible clusters P n . A cluster is a 
partition of (ii, ...in) into Np n ordered subsets p = (ji, ■■■j r ). Then 



>■"■■-" YA\ ( i 1 ( v - 2 °) 



Now from the obvious identity 



we obtain the chain of equations 



6Z 1 S n Z 



$ J nil . . .i n 3Jii---3Ji 1 



(V.21) 



sw 



SJni 1 ...i n Tl. ^ 

We can invert these equations to express the sources as functionals of the connected Green functions, 

and define the master effective action (MEA) as the full Legendre transform of the connected generating 
functional 

Toe [{C r }} = W[{ J n }} - J2 ^ J n }Z II C r ( V ' 23 ) 

n P n p 

The physical case corresponds to the absence of external sources, whereby 

^IWr}] = o (V.24) 
This hierarchy of equations is equivalent to the Dyson- Schwinger series. 

C. The background field method 

The master effective action just introduced becomes more manageable if one applies the background held 
method (BFM) |E3| approach. We first distinguish the mean field and the two point functions 



G\ = <t> 1 (V.25) 



C l 2 J = G» (V.26) 



We then perform the Legendre transform in two steps: first with respect to 4> and G only, and then with 
respect to the rest of the Green functions. The first (partial) Legendre transform yields 

rjA g, {c r }} = r a [0, g, {j n }] - ^ E IT G - ( v - 27 ) 

n>3 ' P„ p 

Here T2 is the two particle-irreducible (2PI) effective action (3l]] 

r 2 [^, G, {J n }} = S[<t>] + ~& k S )jk ~ l - In Det G + J int [<f\ + \& h J intdh + W 2 (V.28) 
and W% is the sum of all 2P1 vacuum bubbles of a theory whose action is 

S'[ip]= % -G- l v 2 + S Q [<p] (V.29) 
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s Q y] = s 1 [4> + ^]- sm - sm 71 ^ - -SjpUj w (v.30) 

where tp is the fluctuation field around <j>, i.e., $ = <fi + (p. Decomposing Sq into source-free and source- 
dependent parts, and Taylor expanding with respect to <p, we may define the background-field dependent 
coupling and sources where 

(T nil ... in = ( m - n )| g "»<i-*»3»+i-jm^' n+1 -^' ro ( V - 31 ) 

m>n 



E (A^...W„ +1 ...i„r'-^ m (V.32) 

m>n 1 



Now, from the properties of the Legendre transformation, we have, for n > 2, 



5J n 1,2 5J n 
Computing this second derivative explicitly, we conclude that 



(V.33) 



dJ„ n\ 2(n-2)! ^ 3 dJ„ dXm 



Comparing this equation with 



we obtain the identity 



sw 2 



Pn P 



£Xnii...»„ n! p _ 



where the * above the sum means that clusters containing one element subsets are deleted. This and 



n>3 P n p r r n>3 P„ p 

allow us to write 

rw^cia}] = + (^'^j - Det G (v- 38 ) 

+tomx»}] -E^*»E II ^ ( v - 39 ) 

n>3 ' P„ p 

This entails an enormous simplification, since it implies that to compute r M it is enough to consider Wi 
as a functional of the Xn, without ever having to decompose these background dependent sources in terms 
of the original external sources. 
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D. Truncation and Slaving: Loop Expansion and Correlation Order 

After obtaining the formal expression for r^, and thereby the formal hierarchy of Dyson - Schwinger 
equations, we should proceed with it much as with the BBGKY hierarchy in statistical mechanics p6[ , 
namely, truncate it and close the lower-order equations by constraining the high order correlation functions 
to be given (time-oriented) functionals of the lower correlations. Truncation proceeds by discarding the 
higher correlation functions and replacing them by given functionals of the lower ones, which represent the 
dynamics in some approximate sense The system which results is an open system and the dynamics 
becomes an effective dynamics. 

It follows from the above that truncations will be generally related to approximation schemes. In field 
theory we have several such schemes available, such as the loop expansion, large N expansions, expansions 
in coupling constants, etc. For dcfinitcncss, we shall study the case of the loop expansion, although similar 
considerations will apply to any of the other schemes. 

Taking then the concrete example of the loop expansion, we observe that the nonlocal \ sources enter 
into W2 in as many nonlinear couplings of the fluctuation field ip. Now, W2 is given by a sum of connected 
vacuum bubbles, and any such graph satisfies the constraints 

nV n = 2i (V.40) 



i-J2 Vn = l ~ 1 ( V ' 41 ) 
where i, I, V n are the number of internal lines, loops, and vertices with n lines, respectively. Therefore, 

Z = 1 + £^K (V.42) 

we conclude that Xn only enters the loop expansion of W2 at order n/2. At any given order I, we are 
effectively setting \ n = 0, n > 21. Since W2 is a function of only X3 to \2i, it follows that the C r 's cannot 
be all independent. Indeed, the equations relating sources to Green functions 



6W1 



have now turned, for n > 21, into the algebraic constraints 

YA\ ( " ( v - 44 ) 

Pn P 

In other words, the constraints which make it possible to invert the transformation from sources to Green 
functions allow us to write the higher Green functions in terms of lower ones. In this way, we see that the 
loop expansion is by itself a truncation in the sense above and hence any finite loop or perturbation theory 
is intrinsically an effective theory. 

Actually, the number of independent Green functions at a given number of loops is even smaller than 21. 
It follows from the above that W2 must be linear on Xn for I + 2 < n < 21. Therefore the corresponding 
derivatives of W% are given functionals of the Xmi m < I + 1. Writing the lower sources in terms of the lower 
order Green functions, again we find a set of constraints on the Green functions, rather than new equations 
defining the relationship of sources to functions. These new constraints take the form 

Xll r = fn(G,C 3 ,...C l+1 ) (V.45) 

p„ p 
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for I + 2 < n < 21. In other words, to a given order I in the loop expansion, only 0, G and C r , 3 < r < I + 1, 
enter into as independent variables. Higher correlations are expressed as functionals of these by virtue 
of the constraints implied by the loop expansion on the functional dependence of W2 on the sources. 

However, these constraints are purely algebraic, and therefore do not define an arrow of time. The dynamics 
of this lower order functions is unitary. Irreversibility appears only when one makes a time-oriented ansatz in 
the form of the higher correlations, such as the 'weakening of correlations' principle invoked in the truncation 
of the BBGKY hierarchy 0. This is done by substituting some of the allowed correlation functions at a given 
number of loops I, by solutions of the /-loop equations of motion. Observe that even if we use exact solutions, 
the end result is an irreversible theory, because the equations themselves arc only an approximation to the 
true Dyson - Schwinger hierarchy. 

To summarize, the truncation of the MEA in a loop expansion scheme proceeds in two stages. First, for 
a given accuracy I, an Z- loop effective action is obtained which depends only on the lowest I + 1 correlation 
functions, say, {0, G, C3, . . . C;+i}. This truncated effective action generates the Woop equations of motion 
for these correlation functions. In the second stage, these equations of motion are solved (with causal 
boundary conditions) for some of the correlation functions, say {Cfc, ...C;+i}, and the result is substituted 
into the I loop effective action. (We say that {Cfc,...C;+i} have been slaved to {4>, G, C3, ...Ck-i}) The 
resulting truncated effective action is generally complex and the mean field equations of motion it generates 
will come out to be dissipative, which indicates that the effective dynamics is stochastic. 



E. Example: the three-loop 2PI EA 



We shall conclude this paper by explicitly computing the 2PI CTP EA for a X(f> 4 self interacting scalar 
field theory, out of the corresponding MEA. We carry out our analysis at three loops order, this being the 
lowest order at which the dynamics of the correlations is nontrivial, in the absence of a symmetry breaking 
background field . 

To this accuracy, we have room for four nonlocal sources besides the mean field and the two point corre- 
lations, namely X3> X4?X5 1 an d Xs- However, the last two enter linearly in the generating functional. Thus 
the three- loop effective action only depends nontrivially on the mean field and the two, three and four point 
correlations. By symmetry, there must be a solution where the mean field and the three point function 
remain identically zero, which we s hall as sume. 

Our first step is to compute Eq. ( V.38), which now reads 



r 4 [G, C 4 ] = ( — )aj (-□ + M 2 ) G ij - - In Det G 



(V.46) 



+ W 2 [(/), {Xn}} ~ ^Xiijkl 



where W2 denotes the sum of 2PI vacuum bubbles of a quantum field theory with quartic self interaction 
and a coupling constant A — xt ( see e Q s - (V.28) and ( V.30D ) up to three loops 



= ( — - ) (A - Xi)aki GVG 



A7 



(V.47) 



+ 



^-XA) m (X-X4) pqr8 G^G kr G ls 



Eq. (VE) yields 



C l i kl = ~i (A - Xi) pqrs G ip G jq G kr G ls 



pqrs 



(V.48) 



Inverting and substituting back in Eq. (V.46), we obtain 
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r 4 [G, <7 4 ] = (-^)cy (-□ + Af 2 ) G y - | In Det G (V.49) 
Xi 3 kiG lJ G kl - f — J XijkiCl jkl 



s~tijkl X s-i—\ ("i—X s~i—\ (~t—V\ /^pqrs 
4 [Lr lp Lx^ l_x fcr U /s J O4 



This functional generates the self consistent, time reversal invariant dynamics of the two and four particle 
Green functions to three loop accuracy. To reduce it further to the dynamics of the two point functions 
alone, we must slave the four point functions. Consider the three loops equation of motion for C4 

[G^G-^G^] CZ S = -ihiki (V.50) 

Solving for this equation with causal boundary conditions yields 

d jm = -i\ pqrs G ip Gi q G kr G ls (V.51) 

(in other words, \i = 0) an d substituting back in Eq. ( V.49| ) we obtain 

T[G] = (^-)cy (-□ + M 2 ) G ij - % - In Det G (V.52) 



Kjkl<- 



f ^ ) \ hk ,G i >-G-"iG l ''G l *\„ 



which is seen to be equivalent to Eq. ( IV.l ). This effective action leads to a dissipative and, as we have 
seen, also stochastic dynamics, which results from the slaving of the four point functions. 



VI. DISCUSSIONS 



In this paper we have introduced a new object, the stochastic propagator G, whose expectation value re- 
produces the usual propagators, but whose fluctuations are designed to account for the quantum fluctuations 
in the binary product of fields. We have introduced the dynamical equation for G which takes the form of an 
explicitly stochastic Schwinger - Dyson equation, and showed that in the kinetic limit, both the fluctuations 
in G become the classical fluctuations in the one particle distribution function, and the dynamic equation 
for G's Wigner transform becomes the Boltzmann - Langevin equation. Each of these results has interest of 
its own. A priori, there is no simple reason why the fluctuations derived from quantum field theory should 
have a physical meaning corresponding to a phenomenological entropy flux and Einstein's relation. 

The notion that Green functions (and indeed, higher correlations as well) may or even ought to be seen 
as possessing fluctuating characters (when placed in the larger context of the whole hierarchy) with clearly 
discernable physical meanings is likely to have an impact on the way we perceive the statistical properties 
of field theory. For example, we are used to fixing the ambiguities of renormalization theory by demanding 
certain Green functions to take on given values under certain conditions (conditions which should resemble the 
physical situation of interest as much as possible, as discussed by O'Connor and Stephens j55|). If the Green 
functions themselves are to be regarded as fluctuating, then the same ought to hold for the renormalized 
coupling constants defined from them, and to the renormalization group (RG) equations describing their 
scale dependence. 

While the application of renormalization group methods to stochastic equations is presented in well-known 
monographs |5q| , our proposal here goes beyond these results in at least two ways. First, in our approach 
the noise is not put in by hand or brought in from outside (e.g., the environment of an open system), as 
in the usual Langevin equation approach, but it follows from the (quantum) dynamics of the system itself. 
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Actually, the possibility of learning about the system from the noise properties (whether it is white or 
coloured, additive or multiplicative, etc.) - unraveling the noise, or treating noise creatively- is a subtext 
in our program. Second, our result suggests that stochasticity may, or should, reside beyond the level of 
equations of motion, and appear at the level of the RG equations, as they describe the running of quantities 
which are themselves fluctuating. 

Indeed, the possibility of a nondeterministic renormalization group flow is even clearer if we think of the 
RG as encoding the process of eliminating irrelevant degrees of freedom from our description of a system 
J57[ . These elimination processes lead as a rule to dissipation and noise, the noise and dissipation in the 
influence action and the CTP-cffcctive action are but a particular case. If the need for such an enlarged 
RG has not been felt so far, the groundbreaking work on the dynamical RG by Ma, Mazenko, Hohenberg, 
Halperin, and many significant others notwithstanding, it is probably due to the fact that the bulk of RG 
research has been focused on equilibrium, stationary properties rather than far- from- equilibrium dynamics 
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VII. APPENDIX: CLOSED TIME PATH CONVENTIONS 

The closed time path (CTP) or Schwinger - Keldysh technique |Q is a bookeeping device to generate 
diagrammatic expansions for true expectation values (as opposed to IN - OUT matrix elements) of certain 
quantum operators. The basic idea is that any expectation value of the form 



(/A |f [4> ..4 {Xn)\ T [4> {x n+1 ) ...<t> m ]\ IN} (mi 



where \IN) is a suitable initial quantum state, x\ to x m are space time points, <p is the field operator, T 
stands for time ordering and T for anti time ordering, may be thought of as a path ordered expectation value 
on a closed time - path ranging from t = — 00 to 00 and back. These path ordered products are generated 
by path integrals of the form 

' DfDcf [0 2 (Xi) ..4 2 (x n ) 1 (x n+1 ) ..4) n ] A s ^y s '^)\ (VII.2) 



where 1 is a field configuration in the forward leg of the path, and 4> 2 likewise on the return leg. These 
configurations match each other on a spacelike surface at the distant future. The boundary conditions at 
the distant past depend on the initial state |/A); for example, if this is a vacuum, then we add a negative 
imaginary part to the mass. We shall not discuss these boundary conditions further, except to note that we 
assume the validity of Wick's theorem (sec [^0[). 

In general we shall use a latin index a, taking values 1 or 2 to denote the CTP branches. Where the 
space time position is not specified, it must be assumed that it has been subsummed within the CTP upper 
index. Also we shall refer to the expression S (4> a ) = S (i/) 1 ) — S* (^> 2 ) as the CTP action. We allways use 
the Einstein sum convention, and if not explicit, integration over space time must be understood as well. 

It is convenient to introduce a CTP metric tensor c a b = diag(l, —1) to keep track of sign inversions. Thus 
CabJ a 4> b — J 1 ^ 1 — J 2 4> 2 . In general, we write an expression like this as J Q a , where J a = c a bJ b ; the index a 
has been lowered by means of the metric tensor. The opposite operation of raising an index is accomplished 

with the inverse metric tensor c ab — (c _1 )° = diag(l, —1). Thus J a = c ab Jb ■ 
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